RENORMALIZED CONTACT POTENTIAL IN TWO DIMENSIONS 



o 

On . R.J. HENDERSON and S. G. RAJEEV 

ON ■ 

Department of Physics and Astronomy, University of Rochester, Rochester, NY 14627 

O 

We obtain for the attractive Dirac 5-function potential in two- 
t^. ■ dimensional quantum mechanics a renormalized formulation that avoids 

reference to a cutoff and running coupling constant. Dimensional 
transmutation is carried out before attempting to solve the system, 
and leads to an interesting eigenvalue problem in iV— 2 degrees of free- 
dom (in the center of momentum frame) when there are iV particles. 
The effective Hamiltonian for N — 2 particles has a nonlocal attrac- 
tive interaction, and the Schrodinger equation becomes an eigenvalue 
problem for the logarithm of this Hamiltonian. The 3-body case is 
examined in detail, and in this variational estimate of the 

ground-state energy is given. 



On 



J|h ' I. Introduction 

> 

Scale invariance, and the ultraviolet divergences for which it is responsible, is 
an essential feature of the quantum field theories, QCD and electroweak, that 
comprise the standard model of elementary particle interactions. The diver- 
gences that are obtained in the process of quantizing the classical field theories 
on which the quantum dynamics are based can be removed via the mathemat- 
ical procedure of renormalization. In condensed matter physics renormaliza- 
tion techniques are used to obtain mathematical models of physical phenomena 
(e.g. phase transitions) which, despite the presence of a characteristic scale, 
provided by a lattice spacing, are scale independent in nature. In contrast, 
renormalization in the context of elementary particle interactions is necessary 
to make the fundamental theory physically sensible. 

That an awkward renormalization procedure is necessary to permit the 
quantum field theory paradigm to successfully describe elementary particle in- 



teractions might be viewed as evidence that quantum field theory is not the 
proper framework for this problem. Some other exotic and finite theory might 
be more conceptually accurate and less mathematically cumbersome. On the 
other hand, one may take Wilson's point of view, [1], on the role of renor- 
malization in particle physics: a renormalizable quantum field theory may be 
viewed as an effective model which approximates for low enough energies (or 
long enough distances) a more fundamental and comprehensive theory. In ad- 
dition, then, to seeking a more fundamental theory underlying renormalizable 
quantum field theories, we may aspire to a deeper understanding of the very 
successful effective theories we have. 

Toward this end we take the point of view that renormalizable interactions 
might be given a finite formulation, [2] , which avoids the need for renormaliza- 
tion altogether, without the necessity of discarding the framework of quantum 
field theory. The renormalized interactions that would be part of this formu- 
lation could not be scale-invariant; the scaling symmetry would be broken ex- 
plicitly at the outset, and not through the renormalization procedure. Other 
properties of the erstwhile "fundamental" interactions might be modified as 
well. 

This approach to finding a finite, effective theory of particle interactions 
would be tantamount to a reordering of the conventional analysis. It would 
require renormalizing the theory completely before embarking on efforts to 
solve it, rather than renormalizing in parallel with the finding of solutions. No 
nonphysical cutoffs or running parameters would appear in the formulation or 
solutions of the renormalized theory. Rather, such a theory could be formulated 
as a well-posed mathematical problem: a set of differential equations with 



appropriate boundary conditions for example. If this point of view is the 
correct one, then the reason we cannot write down such a formulation of QCD 
or electroweak theory is not because they are only low-energy effective theories 
(though they may be) but rather because we have not yet achieved a deep 
enough understanding of these theories to write them down in the simplest 
way. 

Although we believe finding a finite formulation of renormalizable inter- 
actions to be a worthwhile goal, attacking this problem directly, in say QCD, 
appears too formidable at this time. We have in the past, though, found ex- 
amples of simpler, asymptotically free renormalizable theories that lend them- 
selves to this approach. In [3] we renormalized the large- A?" limit of the 1+1 
dimensional non-Abelian Thirring (or Gross-Neveu) model before finding some 
exact solutions. The scale-breaking renormalized interaction manifested itself 
as a restriction on the domain of the Hamiltonian operator. A similar role for 
the Hamiltonian domain was found in [4] in the study of a quantum mechan- 
ical system of two particles in two dimensions attracted by a Dirac 5-function 
potential. The same system was also examined in the path integral picture, 
where it was found that the renormalized interaction appeared as a subtle 
modification to the Wiener measure. 

In this paper we continue with our investigation of quantum mechanical 
particles interacting through an attractive Dirac 5-function potential in two di- 
mensions. We first work out the case of three particles in detail, and eventually 
extend our ideas to the N-body case. 



II. Review of the Two-Body Problem 

Before embarking in the next section on the three-body problem, let us 
recall the renormalized formulation, in the Hamiltonian picture, of the two- 
body problem. In the two-body case, after separating out the center of mass 
coordinate, the original Schrodinger equation in configuration space is: 



where A is the two-dimensional Laplacian, and g is a positive, dimension- 
less coupling constant. We have taken the masses of the two particles to be 
mi = ^2 = 1/2 and chosen units such that h = 1. In momentum space the 
Schrodinger equation reads: 



This eigenvalue problem, however, is nonphysical since, due to scale invari- 
ance, the presence of even a single negative energy solution implies a continuum 
of negative energy states extending down to minus infinity. On the other hand, 
if one attempts to restrict the domain of the Hamiltonian to the positive energy 
sector only, one finds that a complete set of eigenstates cannot be found: there 
will be no zero angular momentum eigenstate, [2]. The Hamiltonian cannot 
be self-adjoint under such a choice of domain. (The situation is analogous to 
the presence of solutions to the Dirac equation of unbounded negative energy, 
and the impossibility of ignoring or "throwing away" these states). 

This system's illness can be cured via renormalization. First one regular- 
izes the system by introducing a momentum cutoff (upper bound). Then the 



2A* A (x) - gS 2 (x)* x (x) = \V x (x) 



(1) 





coupling constant, g, is required to depend on the cutoff in such a way that 
the ground-state energy of the regularized system remains finite as the cutoff 
is removed, i.e. taken to infinity. This procedure removes g from the problem, 
replacing it with a parameter having dimensions of energy, a trade sometimes 
called 'dimensional transmutation'. This new parameter is arbitrary, charac- 
terizes the strength of the interaction, and can be taken to be the ground-state 
energy of the two-body system. 

The customary way to obtain renormalized solutions to the problem, (see 
[5], [6], [7], [8], [9]), then, is to solve the regularized system first, and then 
take the limits of solutions (wavefunctions, scattering amplitudes, etc.) as the 
cutoff is taken to infinity. In [4] we showed that an equivalent, but simpler, 
formulation is given by the following two equations: 

2p 2 V x (p) - lim 2p 2 V x (p) = \V x (p) (3) 

p^oo 

and 

/^( P -)-^ ? )=0 (4) 

where r] X = limp^ 2p 2 ^ x {p). 

These two equations give an example of what we mean by a finite formula- 
tion of a renormalizable theory. The theory at this stage has been renormalized, 
and can be treated as a well-posed mathematical problem. The first equation 
is recognized as a renormalized version of the Schrodinger equation. The in- 
teraction has become the term — = — lim p ^ 00 2p 2 ^(p) on the left hand 
side. This term is the renormalized interaction. Wavefunctions for which r]^ x 



is zero do not take part in the interaction. Interestingly, there is no adjustable 
parameter in the Hamiltonian operator appearing in this equation. In fact, 
the appropriate Hamiltonian in configuration space is the Laplacian with no 
interaction term at all: the deviation from the free theory is contained entirely 
in the specification of the Laplacian's domain. 

The second equation specifies the domain of the Hamiltonian. In con- 
figuration space it is a local condition that implies that wavefunctions with 
zero angular momentum diverge logarithmically at the origin; in configuration 
space, then, the interaction appears as this boundary condition on wavefunc- 
tions. The parameter fx 2 has dimensions of energy, and can be picked arbi- 
trarily. For any choice of (i 2 > 0, the Hamiltonian will be self-adjoint (and in 
fact is a self-adjoint extension of the free Hamiltonian, [10]). Choosing 
corresponds to selecting the strength of the attractive interaction and it turns 
out that —/j 2 is the ground-state energy, and that the free theory corresponds 
to letting [i 2 = 0. 

The above pair of equations can be solved completely, to yield all the 
energy eigenstates in the two-body case, [2] and [4]. Although we cannot 
analytically solve the N-body, or even the three-body, system, we can formulate 
the problem in a similar way, finding the analogs of equations (3) and (4), which 
can then be solved approximately. 



III. Renormalizing the Three-Body Problem 



The original three-body Schrodinger equation in configuration space is, when 
mi = m 2 = TO3 = 1/2 and h = 1: 



(-A 1 -A 2 -A 3 )^ a -^(5 2 (x 1 -x 2 )+5 2 (x 2 -x 3 )+5 2 (x3-x 1 ))*a = A^ A (5) 

where Aj is the two-dimensional Laplacian in coordinate Xj, \I/a = ^a(^i, ^2, ^3), 
and (7 is the positive, dimensionless coupling constant. In momentum space 
this becomes: 

(k* + k% + kl)* x -j^(J d 2 k 12 V x + J d 2 k 23 V x + J d 2 k 3 i^x)=X^x (6) 

where *a = *a(&i, &2, £3) and % = — %)• 

As in the two-body problem, we work in momentum space, in the center of 
momentum frame where K = ki + k 2 + k 3 = 0. Defining pi = §(^(^2 + ^3) — &i) 
(see [11] for a discussion of coordinate choices in the general three-body 
problem), the Schrodinger equation in this frame is: 

(2^3 + \pI)*x ~ J d 2 k 12 V x + J d 2 k 23 ^ x + J d 2 k 31 ^ x ) = X^ x (7) 

where we choose the two independent coordinates to be (k 23 ,pi), such that 
\&A = ^a(^23,Pi), in which case k\ 2 and £31 stand for — §^23 — fpi and — I&23 + 
\p\ respectively. We could just as well have chosen (^31,^2) or (ki 2 ,p 3 ) and 



would have found the same Schrodinger equation, as 2/c| 3 + |p 2 = 2A;| 1 + |p2 = 
2kf 2 + |j>3 since the three particles are identical. From now on we will use this 
symmetry to write this combination simply as 2k 2 + |p 2 . The particle singled 
out by the choice of coordinates (e.g. particle "1" with the choice (A^Pi)) is 
sometimes called the "spectator" particle. Then the other two may be referred 
to as the "interacting pair" , even though in reality all three particles interact 
with each other. 

For the same reasons that (2) was found to be nonphysical, (6) also rep- 
resents an ill-posed physical problem. In this case if one tries to restrict the 
Hamiltonian domain to positive energy states, one will find that states with 
zero angular momentum in the &23 variable are missing, [2] . A renormalization 
of the system is again necessary. By carrying out the renormalization before, 
and independent of, solving the system, we can cast the problem in a form 
analogous to (3) and (4) above. 

The first step in our renormalization program is to regularize the system 
by reconsidering the problem in an artificial momentum space wherein the 
momenta are bounded above by a cutoff. From Wilson's point of view, this 
cutoff represents some scale beyond which we are not entitled to apply our 
low-energy effective theory. Here it is convenient to impose this condition by 
requiring 2k 2 + |p 2 < A 2 . We must also allow the coupling constant, g, to 
depend on A. By doing this we will be able to specify that g change with A 
in such a way that the physics predicted by the effective theory we derive is 
independent of the cutoff. The resulting Schrodinger equation is: 

(2k 2 + \ v 2 - xM = £ J A d% *A (8) 



where v is an arbitrary parameter with dimensions of momentum, J indicates 
the the restriction 2k 2 + |p 2 < A 2 , the sum has three terms, and from now 
on fei, &2 and ks will correspond to, and be used interchangeably with, k 2 3, &3i 
and ki2 respectively. 

The presence of A breaks scale invariance, such that there is no longer the 
instability problem preventing us from allowing negative energy states. Also, 
if we define f^iipi) = ^j£^w f A d 2 ki^^, then (8) can be written: 

(2k 2 + 3 -p 2 -\M = J2fx,^) ( 9 ) 

i 

If we decide to consider to be a function of fc 2 3 and p\, then it is also 
useful to note: 

P2=^(h + k 1 )-k 2 ) = -k 2 3-^Pi (10) 



and 



21-- - - 1 _ 



(11) 



It is left to determine what g(A/u) must be in order that this regularized 
system will have a well-defined limit as we let A — > oo. By integrating (8) over 
any of the ki it is not difficult to show, [2], that g(A/u) must satisfy, for large 
enough Ajv: 



SW") Jk<A 2fc 2 + f 2 

where we keep in mind that v 2 is arbitrary. (12) implies that g(A/u) — > as 
we let A — > oo, indicating that the interaction we consider is asymptotically 
free. This expression for g is, with v 2 replaced by fx 2 , precisely what was found 
during the renormalization of the two-body system, [4]. This is a result we 
should expect. In fact it would be have been reasonable to fix g(A/v) using 
the two-body result rather than to rederive it in the three-body case. This 
correspondence allows us to identify — v 2 = —fi 2 as the binding energy of the 
two particle bound state, i.e. the ground-state energy of the two-body system. 
From now on we will use —fx 2 to denote this quantity. 

By taking the /C23 — > 00 limit of (8) we can make the identification: 



lim 2k 2 ^{ = 9 -j^ [ A rfWJ = fUPi) (13) 



^23 



This equation has counterparts with /c 2 3 replaced by £31 or k±2, such that (8) 
could also be written: 



(2^ 2 + |p 2 - AM = V lim 2kM (14) 



This equation does not contain g, and will remain valid as we let A — > 00. 
In this limit, (14) becomes the three-body renormalized Schrodinger equation, 
comparable to (3) in the two-body case. Again, no coupling parameter appears 
in the renormalized Hamiltonian. The parameter characterizing the interaction 



strength, which we continue to take to be the two-body binding energy, will 
appear in the domain (boundary condition) equation. 

We can find the domain equation by combining (12) and (13). Together, 
as A — > oo, these equations impose the following conditions on wavefunctions: 



where f\(pi) = lim k .^ 00 2kf^> \. One way of understanding the choice of 
g(A/v) = g(A//j) in (12) is that g must have the form which ensures that 
(13) will continue to hold true as we let A — > oo. 

(15) is in fact three equations, one each for i = 1, 2, 3. They are integrals 
over momentum space, but in configuration space they are local conditions 
dictating that the wavefunctions either go to zero or diverge logarithmically, 
at a rate determined by fi, whenever the coordinates of two particles are made 
to coincide. That the wavefunction can blow up under these circumstances is 
unusual, but not proscribed since the singularities are square-integrable. 

We thus find that the three-body analogs of the two-body equations (4) 
and (3) are (15) and the three-body renormalized Schrodinger equation: 



The equations (15) and (16) comprise the renormalized, finite formulation 
of the three-body problem we seek. The form of these equations indicates 
that the unknown parts of the wavefunctions are essentially their limits as 
the relative momenta of particles are taken to infinity, i.e. the functions f\^, 




(15) 




(16) 



i = 1,2,3. We will require the particles to satisfy Bose statistics, in which 
case must be symmetric under permutations of particle indices. Then the 
fx,i must all be the same function, which we will call fx- fx is a function of 
one momentum variable rather than two, and so the eigenvalue problem has 
already been simplified. The analogous quantity, ry«i< A = \iva p ^ 00 2p 2 '^> x(p), in 
the two-body case was simply a constant which could be fixed by normalization, 
so that the eigenvalue problem could be trivially solved. In the N-body case 
we will find an eigenvalue problem in N — 2 degrees of freedom in the center 
of momentum frame. 

Although the crux of our problem in solving the system is finding the 
functions fx for all eigenvalues A, (15) and (16) contain not only fx but the 
wavefunction ^x as well. It would be convenient to eliminate ^>x from the 
problem altogether, and obtain the equation satisfied by fx alone. This we can 
do by employing the Lippmann-Schwinger formulation, [12] , of the Schrodinger 
equation. The Lippmann-Schwinger approach is conventionally reserved for 
scattering states, but here we will find it a convenient starting point in both 
the scattering and bound state sectors of the theory. 

The Lippmann-Schwinger equation for our three-body system can be writ- 
ten down by inspection using (16). It is: 



We choose for concreteness to write (17) in the basis (£23, pi), g (pi, fas) is 
arbitrary, excepting the fact that it must be symmetric under permutations of 
particle indices. 



^x = g(pi,k 23 )5(2k 2 + ^p 2 -X) + 



"Lift (Pi) 



(17) 



2k 2 + \p 2 - A=Fze 



The energy e > in (17) is infinitesimally small. It must appear when 
A > to make the division by the singular operator 2k 2 + |p 2 — A well- 
defined. The choice of sign, =F,m the denominator corresponds to a choice 
of boundary conditions: the upper sign implying outgoing scattered waves 
(which is ordinarily the physical case) , and the lower sign meaning converging 
scattered waves, e can eventually be taken to zero, but will serve to regulate 
otherwise divergent integrals in intermediate steps. We add the label ± to 
to signify this choice of boundary conditions. 

The first term in (17) represents the kernel of the operator 2k 2 + \p 2 — 
A, which appears on the left hand side of (16). It is a solution to the free 
Schrodinger equation. In the case A > it represents the unscattered portion 
of the wave, and in configuration space gives the wavefunction its asymptotic 
behavior as particles become infinitely separated. When A < this term 
disappears. Negative energy states are comprised of only the second term, and 
their form is specified by (17) up to the unknown function f\. All information 
contained in (16) is also in (17). 

At this point we must part ways with the usual Lippmann-Schwinger 
analysis. No potential appears in (17), and we cannot use the self-consistency 
of this equation alone to solve for the wavefunction. In the place of usual term 
involving the potential stands the unknown function f^. The only further 
information we can glean about from (17) comes from taking the limit of 
this equation as /C23 — > oo. Since f^ipi) = h m A;23^oo ^k^^x, taking this limit 
and using (10) and (11) gives: 



lim (fx(k23) + fx(-h3)) = 



(18) 



However, the bosonic symmetry requires that /a (—^23) = fx (^23), so that 
(18) becomes the simple asymptotic condition: 



lim f x (k 23 ) = (19) 

fc 2 3^0O 



To gain more information about fx, we need to use (15). This once more 
highlights the fact that in this system, the interaction between particles is 
encoded in the boundary conditions on wavefunctions rather than in a conven- 
tional potential energy term. Inserting (17) into (15) and performing integra- 
tions where possible yields the following implicit equation for fy. 



(W \ \\--p 2 \)^e(\--P 2 ))ft(p)-- fd 2 k — 

= ^Jd 2 k g(p,k)5(2k 2 -^p 2 -X) 



(20) 



This ungainly looking equation contains all the information about wave- 
functions and energy eigenvalues we need. Imposing the requirement that f^ 
has asymptotic behavior such that the integral in (20) converges will ensure 
(18) holds true. When A > 0, this equation is an inhomogeneous linear equa- 
tion in f^. Operating on on the left hand side of (20) is a symmetric 
integral operator, which must be inverted to get for any given g. When 
A < 0, the inhomogeneous term proportional to g vanishes identically, leaving 
an eigenvalue problem for f^ and A. In the following section we partially solve 
the problem of finding f^ for positive and negative energies, and show that 



the negative energy sector yields an explicit form for a renormalized, nonlocal 
effective Hamiltonian. 



IV. Toward Solving the Three-Body Problem 



All solutions to the free Schrodinger equation have positive energy and can be 
categorized by the following choice of basis: 

^ ree = g(p 1 ,k23)S(2k 2 + ^p 2 - X) (21) 

where 

g{piMz) = ^2s 2 (k 2 3 - ^(«n(2) - ^n(3)))£ 2 (Pi - («n(2) + «n(3))) (22) 
n 

The constant vectors k±, R2 and K3 are the momenta of the three particles 
and label the states. There are no restrictions on them other than they repre- 
sent the particles in the center of momentum frame wherein R\ + R2 + R3 = 0. 
The energy is A = k\ + k 2 , + k\. The sum in (22) is over the six permutations 
of (1, 2, 3), and gives g the bosonic symmetry we require. For our purposes it 
is convenient to absorb the /C23 dependence into pi using the delta function in 
(21) and to diagonalize the angular momentum in the angle, 6*23, of ^23- In 
this case g becomes: 

£n(Pl,023) = 

e in9 23 Y^ e -ine m(3)5{k 2^ _ 1 (-^ _ Rm) f )S 2 (p 1 - (« n(2 ) + «n(3))) 

n 

(23) 



where n is an integer, and $ij denotes the angle of — Rj). 

Recall that in the two-body problem only those states carrying zero angu- 
lar momentum participated in the interaction. The wavefunctions with nonzero 
angular momentum were just the free ones. We find a similar situation in the 
three-body case. The free wavefunctions 

*n,A = 9n(Pu0 23 )d(2k 2 + ^p 2 - A) (24) 



for n ^ are the positive energy solutions we seek in the sector of the Hilbert 
space where the angular momenta of any pair of particles is nonzero. Clearly 
(24) solves (16) with f\ = 0, for any n. For n = 0, due to the integration 
over 6>23, these wavefunctions also satisfy (15) for f\ = 0. Note that the right 
hand side of (20) becomes zero in this case. This is consistent with the linear 
operator on the left hand being nonsingular. 

Free wavefunctions with quantum number n not equal to zero do not 
scatter. Nontrivial, i.e. interacting, wavefunctions can therefore be taken to 
have the form (17), with g(pi, ^23) — ► 9oi.Pi)- 9o is the symmetrized (goipi) = 
9o{p2) = 9o{Pz)) function, for n = 0, given in (23). It is independent of 6*23, so 
the integral on the right hand side of (20) is trivial, and this equation becomes: 



= -e{\-\p 2 )g Q {p) 

7T A 



(25) 



Positive energy scattering solutions, exist for all A > and all go. Note 



that the operator acting on in (25) is rotationally invariant. One can let 
fx ip) ~ e%l&v > where the integer / is the total angular momentum of the three 
particle system, and solve each angular momentum sector independently. It 
is also worth noting that a rotation in the angle 6*23 is not a symmetry of the 
interacting Hamiltonian; although scattered wavefunctions have a free part 
independent of this angle, the scattered components of these states will in 
general depend on #23. Equation (25) is as far as we will take our analysis in 
the positive energy sector. 

For the consideration of negative energy states, let us take A = — if. In 
this case, (17) gives the form of the wavefunctions as: 



= 2k* + P + if (26) 



The infinitesimal parameter e can be set to zero in this case. When A = — rj 2 
(25) becomes an eigenvalue problem for f\ and A: 

K7W + lp 2 ))h(P) -If d 2 k 2 Mk) = (27) 

V ^ 2 t TV J (fc2 +p 2 +p . k+ !L) 

For a separable attractive potential (i.e. one having the form < k\V \k' >= 
—v(k)v(k')) the Faddeev integral equation for the T-matrix, [13], written as 
an equation for the wavefunction, can be reduced to an equation in just one 
variable. Bruch and Tjon, [14], have shown that, in the limit in which such 
a potential is made to have zero range in real space, this equation is precisely 
the eigenvalue equation we have found in (27). 



Our eigenvalue problem can be brought to a more conventional form by 
writing it in terms of dimensionless variables. We let p = rjx, k = rjy, and 
fx(p) = u\(x). Then (27) becomes 

ln(l + f * 2 W*) --! <?y U 2 f} t = H^)u x (x) (28) 

V 2 / n J (y z + x 2 + x ■ y + ±) r] z 

Denoting the operator on the left hand side of (28) by W, we can rewrite 
that equation as: 

Wu x = ln(4)«A (29) 
We find, therefore, a linear eigenvalue problem for the eigenfunction u\ 

2 

and the eigenvalue ln(^). Note that W is a symmetric, nonlocal operator. Its 
eigenfunctions, u\, give us the undetermined part (see (26)) of the wavefunc- 
tion, * A . 

Since the energy is — rf , can identify the operator 

H = -^e- w (30) 

as the Hamiltonian, with one degree of freedom (in addition to the center of 
mass) effectively integrated out, of the renormalized system in the negative 
energy sector. We thus find that the logarithm of the effective Hamiltonian of 
the spectator particle, in momentum space, is an integral operator. The form 
of the integral operator W reveals that the renormalized, effective "spectator 
Hamiltonian" includes an attractive potential which is nonlocal. It is important 



to realize, however, that the nonlocality of this Hamiltonian arises because we 
have effectively integrated out one degree of freedom, a simplification we could 
make because the actual interaction is of zero range. A positive attribute of 
the formulation in terms of this one-particle Hamiltonian is that the boundary 
conditions, and therefore the interactions, are built in, and do not come in 
through some supplementary condition on wavefunctions. 

The operator W consists of two terms: an unbounded positive kinetic- 
energy-like operator (a multiplication by ln(|x 2 + 1) in momentum space) and 
an interaction term (the integral operator in (28)). In light of numerical and 
variational evidence we conjecture, but have not proved, that the interaction 
part of W is bounded below. Certainly the Hamiltonian including both the 
"kinetic" and "potential" energy terms is bounded below. For a proof of this 
see [14]. 

If the eefective potential energy term alone is in fact bounded below it 
indicates that there are a finite number of normalizable bound states, as well as 
a continuum of negative energy scattering states wherein a bound state of two 
particles scatters from the third particle, and the kinetic energy of scattering is 
less in magnitude than the two-particle binding energy, — . These states will 
have counterparts in the positive energy sector for which the kinetic energy 
dominates. Also in the positive energy sector will be scattering states of three 
unbound fundamental particles. In fact Bruch and Tjon, [14], have carried out 
a numerical study of (28) and concluded that there are just two three-particle 
bound states, a result consistent with the above picture. 

The expression for H in (30), in combination with the explicit form of W, 
is one of our key results. It gives us a rare glimpse at an explicit renormalized 



Hamiltonian, that in most theories can at best be asymptotically approached 
order by order in perturbation theory. In this system the entire interaction 
appears in the domain of the renormalized Hamiltonian. The fact that the 
interaction is of zero range allows us to obtain a compact form for the eigenvalue 
problem by effectively integrating out one degree of freedom. That the simplest 
way to write the resulting Hamiltonian is in terms of its logarithm is interesting, 
but should not surprise us in light of the prevalence of logarithmic dependences 
in asymptotically free theories. We will find in the N-body case as well that 
the Hamiltonian is best written in terms of its logarithm. Before doing this, 
however, we will analyze the three-body case a little further, finding in the 
next section an approximation to the ground state and ground state energy. 



V. Approximating the Three-Body Ground State 



Our system of three particles is invariant under rotations, so that we may use 
this symmetry to simultaneously diagonalize the Hamiltonian and the total 
angular momentum. Letting the integer / be the angular momentum quantum 
number, the form of negative energy wavefunctions becomes: 

* X > 1 ~ 2k* + Ip* + v 2 ( } 

where 9i is the angle of pi. 

Letting u\j(x) = f\j(p), the eigenvalue equation, (28), can be decom- 
posed into separate angular momentum sectors, giving us: 



3 f 00 dy U(x + y + ^-xy-(x + y+l)) 

• /() { - >!j) sj(x + y + \) 2 - xy 



H 2 



(32) 



The ground state will correspond to the (or one of the) / = negative 
energy state(s). Letting — 7 2 be the ground state energy, the form of this 
state, which we will call \l/o 5 is: 



,T/ — Zli fo(Pi) /oqN 

*° " 2FWTf (33) 



with fo(p) = uq(x) being the eigenstate corresponding to the lowest eigenvalue 

2 

ln(^a-) solving the equation: 



ln(|s + l)u(x) - 2 r dy^ ^ = H^)u{x) (34) 

By definition u solves this equation for the value rf = 7 2 . Equivalently we 
can regard uq as the function which minimizes the quadratic functional: 



Q[u] = [°° ln(L + l)u 2 (x) -2 dx dy < x )<y ) (35) 

^0 2 ^0 io J(x + y + 1)2 _ xy 

subject to the constraint J °° dxu 2 (x) = 1. 

One way of obtaining an approximation to the ground state wavefunction 
is to numerically solve a discretized version of (34). That is, we can diago- 
nalize the matrix version of the operator W, Wij, obtained by replacing the 
continuous variables < x, y < oo by discrete ones < iA, jA < NA for 
integers i,j and iV and positive A. The task then is to diagonalize Wij. Doing 
this amounts to making a discrete approximation to the regularized problem 
within which there is a high momentum cutoff. For large iVA the unique 
positive eigenvector should give a useful picture of the ground state of the 
renormalized problem. 

Using this method to obtain a numerical estimate of the ground state 
wavefunction shows that uq has a long, power law type tail. In fact uq(x) is 
quite well approximated by a function of the form y/b/(b + x). This rational 
function is normalized such that the integral of its square on the positive real 
line is unity, and therefore can be taken as a variational ansatz for uq. Inserting 



this form in the quadratic functional Q[u] and minimizing with respect to the 
parameter b yields the upper bound on the ground state energy: 

A = - 7 2 < - /U 2 e 2 - 6 « -13.5// 2 (36) 

Recalling that — /j 2 has the meaning of the two-body ground state energy, 
we see that the three-body energy is significantly less, which is physically 
reasonable. Also, this bound is consistent with the estimate Ao ~ — 16.1/x 2 
found by Bruch and Tjon using a numerical diagonlaization of the Hamiltonian. 
In addition they find one other bound state, with spectator particle angular 
momentum / = 1, having energy m — 1.25/i 2 . 



VI. The N-Body Problem 



The original Schrodinger equation of the N-body problem in momentum space 
is: 

N r 

£ d% ^ x = a * a (37) 

i=l ^ ' i<j J 

where we again take fi = 1 and the particle masses to be m; = 1/2. kij = 
\(Jz% — kj) are, as in the three-body case, the relative momenta of pairs of 
particles. As before the system is ill-defined due to scale invariance. We 
regularize the problem with a momentum cutoff A such that £V k 2 < A 2 , and 
define for each pair of particles the functions: 

ft* = < 38 > 

g(A/n) is as given in (12), and the superscript A on the integral indicates the 
restriction £V k\ < A 2 . Recall that —fx 2 is the two-body ground state energy. 

As A is taken to oo in (38), g{K/ fj) is driven to zero, and picks out only 
the logarithmically divergent part of the integral over the relative momenta of 
particles. The renormalized Schrodinger equation then becomes: 

JV 

(X;*f-A)* A =x;/ Alii (39) 

i=l i<j 

The functions f\^j are the limits as A — > oo of the f\ij, and in this 
limit (38) and (12) ensure that f\^j = limfe i;) ._ >00 2k 2 ^\. f\^j is independent, 
then, of kij, and if we continue to consider the particles to be bosons, then the 



necessary permutation symmetry will imply that only one of these functions 
need be specified. The others will follow by permutation of particle indices. 

In a similar manner as in the two- and three-body problems, we can show 
that the renormalized Schrodinger equation, (39), must be supplemented by 
the domain equations: 



/ 



^(*A"^f^) = (40) 



There is no difficulty in using the Lippmann-Schwinger form of (39) in 
combination with (40) to derive an equation like (20) for the N-body problem. 
The full equation is not very illuminating, however, so we do not reproduce it 
here. Rather we confine ourselves to the negative energy sector and again let 
A = — i] 2 . The wavefunctions then have the form: 



Ef=i + v 2 



*A = JT! A '\ (41) 



Inserting this form into (40), taking, for example, the integration to be 
over the relative momentum ki2, and transforming to dimensionless variables 
Xi = ki/r], yields the eigenvalue equation: 



ln(l + 1^+ i2 f + £ - I / d Wffi' 1 * = ln(fj)«,« 

i=3 J \2^i=l x i + *■) ' 

(42) 



where u\ j ij(xi, xn) = f\,ij(ki,---,kN)- One must keep in mind that u\^j 
is independent of \{xi — Xj), so that in the center of momentum frame this 



function depends on N — 2 two-dimensional vectors. The linear operator on the 
left hand side of (42) is the N-body version of the operator W. We may, as in 
the three-body system, identify the renormalized spectator Hamiltonian in the 
negative energy sector as the exponential of this operator, i.e. H = —fj, 2 e~ w . 

We conclude that the N-body problem can be given a renormalized for- 
mulation that is in principle no more complicated than for the three-body 
problem. Of course it is much more difficult to solve, or even to estimate its 
ground state and ground state energy. Our most intriguing result, however, 
that the simplest renormalized formulation is in terms of the logarithm of the 
Hamiltonian which is an integral operator in momentum space, continues to 
hold in the negative energy sector of the N-body problem for arbitrary N. 

We have found, therefore, for all N, a renormalized, i.e. finite, formulation 
of the system of N particles interacting through an attractive Dirac 5-function 
potential. 

An undesirable feature of our formulation is the schism between the treat- 
ment of positive and negative energy states. States for all positive values of 
energy exist, and are found by inverting a linear integral operator in momentum 
space. Negative energy states, in contrast, come as solutions to an eigenvalue 
problem. This is the price we pay for combining equations (16) and (15) into a 
formulation of the problem in which the boundary condition on wavefunctions 
is built in, and not a supplementary condition. 

In a sense, however, all the essential features of the system are included 
in the negative energy sector. In this sector are contained normalizable bound 
states, scattering amongst these composite particles, as well as scattering of 
the composite particles with the elementary excitations of the theory. It is 



possible, in fact, that as we allow N — > oo, keeping the ground state energy 
finite, that it is only this sector which evolves into the resulting field theory. 
The field theory found this way should admit a finite formulation along the 
lines presented here, and is a natural direction for future research. 



Acknowledgments 



This work was supported in part by funds provided by the U.S. Department 
of Energy under grant DE-FG02-91ER40685. 

References 

1. K.G. Wilson, Phys. Rev. D2, 1438 (1970); Phys. Rev. B4, 3174 and 3184 
(1971); Phys. Rep. 12C, 75 (1974); Rev. Mod. Phys. 55, 583 (1983). 

2. R.J. Henderson, Ph.D. Dissertation, University of Rochester, Rochester, 
NY (1997). 

3. R.J. Henderson and S.G. Rajeev, Intl. J. Mod. Phys. A10, 3765 (1995). 

4. R.J. Henderson and S.G. Rajeev, to appear in J. Math. Phys. (1997). 

5. C. Thorn, Phys. Rev. D19, 639 (1979). 

6. K. Huang, Quarks, Leptons and Gauge Fields, World Scientific, Singapore 
(1982). 

7. K.S. Gupta and S.G. Rajeev, Phys. Rev. D48, 5940 (1993). 

8. C. Manuel and R. Tarrach, Phys. Lett. B328, 113 (1994). 

9. G. Amelino-Camelia, Phys. Lett. B326, 282 (1994); Phys. Rev. D51, 2000 
(1995). 

13. R. D. Amato and J. V. Noble, Phys. Rev D5, 1992 (1972). 

14. L. W. Bruch and J. A. Tjon, Phys. Rev. A19, 425 (1979). 

10. S. Albeverio, F. Gesztesy, R. Hoegh-Krohn, and H. Holden, Solvable Models 
in Quantum Mechanics, Springer- Verlag, New York, NY (1988). 



